Dynamical excitonic effects in metals and semiconductors 
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The dynamics of an electron-hole pair induced by the time-dependent screened Coulomb interac- 
tion is discussed. In contrast to the case where the static electron-hole interaction is considered we 
demonstrate the occurrence of important dynamical excitonic effects in the solution of the Bethe- 
Salpeter equation. This is illustrated in the calculated absorption spectra of noble metals (copper 
and silver) and silicon. Dynamical corrections strongly affect the spectra, partially canceling dy- 
namical self-energy effects and leading to good agreement with experiment. 
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The dynamics of an electron-hole pair in a many body 
system is a strategic field of research with applications 
to many different aspects of solid state physics 0. For 
a large number of systems the observed light-absorption 
spectra largely deviate from independent-particle calcu- 
lations 23. When carried out including self-energy ef- 
fects, spectra calculated for semiconductors and insula- 
tors are characterized by a shift toward high energies with 
respect to experiment. Those deviations are corrected by 
including the electron-hole interaction 0]. The strength 
of these modifications increases as the inverse of the di- 
electric constant of the system. In insulators, where the 
electron-hole interaction is only weakly screened, sharp 
peaks with energy below the optical gap (bound exci- 
tons) can be observed in the experimental spectra 3]. If 
this description is extrapolated to the metallic case the 
natural conclusion is that the electron-hole interaction 
has a negligible effect on the optical spectra of metals as 
the static electron-hole interaction is completely screened 
by the long-range part of the dielectric function. This 
simple argument has been considered definitive to assert 
that there are no excitonic effects in metals Q. In sil- 
ver and copper, however, independent -particle calcula- 
tions overestimate the experimental absorption spectra 
strengths by ~ 30% |E 13 1 a deviation that may be 
explained only in terms of many-body effects beyond the 
independent-particle approximation. 

The standard approach to account for the electron- 
hole interaction in optical spectra is based on the so- 
lution of the Bethe-Salpeter equation (BSE) P, ||| for 
the two-particle Green's function. An important ingre- 
dient of the BSE is the electron-hole interaction, de- 
scribed by the screened, time-dependent Coulomb inter- 
action W {Y^r'\t — t'). However, the BSE with a time- 
dependent interaction is hardly solvable and, as com- 
mon practice, the electron-hole interaction is assumed 
to be instantaneous; this is equivalent to approximat- 
ing the time Fourier transform of W with its static 
value, W {yi,Y2tOJ = 0). This approximation is verified 



a posteriori through the comparison with the experi- 
ment and physically corresponds to the assumption that 
the electron-hole scattering time is much longer than 
the characteristic screening time of the system (roughly 
speaking, the inverse of the plasma frequency). Indeed, 
the static approximation is expected to work well for 
transition energies much smaller than the plasma fre- 
quency 0. However the most striking examples of sys- 
tems that do not fulfill this condition are silver and cop- 
per. The well-known sharp plasmon of silver, which dom- 
inates the electron-energy-loss spectra (EELS) 0,0) is 
located just above the interband gap (^3.9 eV). Similarly, 
the EELS of copper shows strong, broad peaks in the op- 
tical range 1^. 

In this letter we prove that dynamical excitonic effects 
in copper and silver are possible when the BSE is solved 
with a time-dependent electron-hole interaction. Their 
experimental optical spectra are correctly explained only 
thanks to a delicate interplay between dynamical exci- 
tonic and self-energy effects. In order to clarify the re- 
lation between the present approach and the previous 
results obtained within the static approximation 0], we 
show that the optical spectrum of silicon can be obtained 
without using the approximations commonly employed in 
the static approach, that is without neglecting quasipar- 
ticle renormalization factors and dynamical electron-hole 
interaction effects. 

The absorption spectrum is given by the imaginary 
part of the dielectric function e (ut) = 1 — SttA'I'P (oj) A, 
where P (w) is the matrix representation of the polar- 
ization function in the non-interacting electron-hole ba- 
sis and A is a vector embodying the corresponding op- 
tical oscillators. The single-particle states are calcu- 
lated by means of density-functional-theory (DFT) in 
the local-density-approximation (LDA) [loj , while quasi- 
particle (QP) corrections are added on top of the DFT- 
LDA band structure following the implementation of the 
GW method ^ described in Ref. The polariza- 

tion function is obtained by solving the BSE, an integral 
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equation for the four point electron-hole Green's function 
L (ti, ^3, ^4)- As we are interested in the polarization 
function P (t) = —iL {t, 0; t, 0), the BSE can be rewritten 
as [13: 

P (<) = P(°) (t) - J dh P(") {t - ii) VP (<i) 



J I dtidt2 L(°) (i, t2] t,ti)W (ti - t2) L (ti 



0;i2,0), 



(1) 



where the BSE kernel is decomposed into a sum of the 
instantaneous bare electron-hole and exchange interac- 
tion V and of the time dependent screening contribution 
W: W{t) = V5{t) +W{t). Here we have explicitly 
quoted only time variables, since time dependence is the 
main concern of this work. L^*^) is the non-interacting 
electron-hole Green's function whose matrix elements are 

^(ci)k,c'i>'k') (^Ij ''Z; ^3, ^4) = Sy^y'6c,c'Sk^k' Zc]iZyk 
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c, V, k being conduction, valence band and k-point in- 
dexes, ii'nk (with n = c,v) and Znk (smaller than 1) 
are the QP energies and renormalization factors, respec- 
tively. The latter represent the weights of the QP peak 
in the many-body single-particle spectral function. The 
more Z„k differs from 1, the more the high energy struc- 
tures in the spectral function (like plasmonic replicas) 
become important, due to the coupling of the QP with 
system excitations. Those high-energy peaks are not vis- 
ible in the optical energy range but, nevertheless, sub- 
tract intensity from the QP peaks. Very little is known 
about the role played by the Z factors in optical spec- 
tra calculations [lj| . If the Z factors are included in an 
independent-QP calculation or even in the BSE (see 
below) , the intensity of the resulting spectra is strongly 
underestimated, both in metals and in semiconductors. 
Thus the Z factors are commonly set to 1 by hand in 
the solution of the BSE or in the calculation of the 
independent-QP spectra. This is, again, an approxima- 
tion, needed to reproduce the experimental results, with- 
out a sound theoretical justification. 

Because of the time dependent term W(ti— ^2), 
Eq. 1^ cannot be rewritten in terms of P (t) only. For 
this reason the BSE is considered hardly solvable (if not 
"practically unsolvable" (14| ) and for computational con- 
venience the electron-hole interaction is approximated as 
static, W (t) w W (tj = 0) S{t) 0, 0. Thus Eq. ^ can 
be formally solved by means of a Fourier transform: 



P (lu) = P(") (w) - P(°) {uj) f V + w) P [lu 



(3) 



This is the static BSE (SBSE) commonly applied neglect- 
ing the renormalization factors in Eq. |(2Jl, i.e. taking 
Znk = 1- It yields optical spectra in good agreement 



with experiments in semiconductors and insulators 
When applied to copper and silver however, the SBSE 
result (dotted lines in Figl^J is indistinguishable from 
the independent-QP calculation, without improving the 
agreement with experiment. The inclusion of the appro- 
priate Z factors in the independent-QP polarization, P°, 
leads to strongly underestimated absorption spectra in 
all cases Hence, we look for a solution of equation 
(1) without the two major approximations employed in 
the SBSE, i.e. keeping the Z's smaller than 1 and W 
frequency dependent. To this end we expand L'^'WL in 
powers of W. Using generalized indexes K := (cuk) the 
first order term of this expansion, P*^^^ (t), is given by: 



P^L it) 
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From Eq. |2l it is straightforward to see that 

l'^I (i, t2;t,ti)=i [P^l (t-ti) e'^"!--! (ii - t2) 

+^Ki (^-i2)e-*^=i'^i(*^-*^)0(t2-ii)] , (5) 

(ii,0;t2,0) [4°) {t2)e-'^'^^-^'^'^-*-'>0{t,~t2) 
+P^l{ti)e'''^2-2^*'-''^e{t2-t,)\ , (6) 

that inserted in Eq. Q), casts P*-^-* (t) as a time convolu- 
tion of three terms (as shown diagrammatically in Fig^ 
left diagram). 
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FIG. 1: Diagrammatic representation of the first (P^^^ (t) for 
ti > t2, left diagram) and second (p(-^'") (t), right diagram) or- 
der contributions to the polarization function P (t) according to 
the BSE (see text). Crosses indicate the time points where the in- 
coming and outgoing non— interacting Green's functions are "cut" 
according to Eqs. 15161 . 



As a consequence, in the frequency domain P*^^) {uj) 
has the form: 



p(i) (w) = -p(0) (u;) lu^"^ (uj) + nC*) (iu)] P(") (u;) , 



(7) 
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ScikJ, W(+) (u;) be- 



ing the Laplace transform of W(i). The two terms 
denoted by (a) and (5) correspond to the two possi- 
ble time orderings of the interaction ends (ti > <2 for 
term (a), shown in Fig^ t2 > ti for term (b), not 
shown). Eq. {Tj) can be thought of as the first order ex- 
pansion of P (a;) in the frequency-dependent interaction 
n (w) ^ n^"") (w) + n('') (w), which replaces W of the 
SBSE. Thus a partial summation of the BSE can be per- 
formed writing: 

P (w) = P(°) {u) - P(°) (lu) [V + n (lu)] P (u;) . (8) 

This is the Dynamical Bethe-Salpeter equation (DBSE), 
the central result of this work. The diagrams summed 
up in Eq. JHJ are those containing the ladder series of re- 
peated electron-hole interactions with non overlapping 
(in time) interaction lines. The poles of P {uj), 0a, will 
be given by the solution of the equation [p(°) (0a)] + 
V -I- n (0a) — 0. In contrast to the kernel of the 
SBSE, n(0A) is not hermitian and, consequently, 0a 
is in general complex. Its imaginary part gives the in- 
verse excitonic lifetime. Thus the interacting electron- 
hole states are actually dressed excitons, or quasiexci- 
tons. This agrees with what has been already found 
in the core exciton limit and emphasizes the anal- 
ogy between the DBSE and the Dyson equation. Con- 
sequently, as in the single-particle problem, we expect 
to find similar renormalization effects on the quasiexci- 
tonic Green's function. To develop further this aspect we 
expand linearly the smooth function W'-^^ {lu) around 
the non-interacting electron-hole energies, obtaining 
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0KiK2 (^ — £'C2k2 + Ey^\^^). 

is the static limit 
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of the dynamical Bethe-Salpeter kernel which turns 
out to be quite similar to the kernel of the SBSE. 
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are the exci- 



tonic dynamical-renormalization factors. Thus Eq. (jSJ 
can be strongly simplified in the case of noble metals 
where the effect of n*^**^ -I- V is very small. The corre- 
sponding polarization function P (cj) is approximatively 
given by: 



-1 

K1K2 



uj-E, 



C2^2 
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(9) 



with — ^ciki^'uiki<^KiK2- The connection be- 

tween dynamical excitonic and self-energy effects is now 
clear. = 1 — /3„k, where the negative factor /3„k, 

the frequency derivative of the self-energy, is the weight 
lost by the QP because of the coupling with the excita- 
tions of W {ijj). The excitonic factors 0, instead, are due 
to the modification of such coupling as a consequence of 



the electron-hole interaction. Those two effects tend to 
cancel each other but the cancellation is, in general, not 
complete, as exemplified in Fig|21 for copper and silver. 
The SBSE calculation (dotted line), with Znk = 1 and 
= 0, overestimates the experimental intensity (cir- 
cles), while the inclusion of the Z's only (dashed line) 
underestimates it 0. In the DBSE (full line), obtained 
solving Eq. 0, the dynamical factors partially com- 
pensate for the Z^'' factors yielding a spectral intensity 
in good agreement with experiment. 
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FIG. 2: Absorption spectrum of bulk copper, silver, and silicon. 
Dotted line: SBSE without renormalization factors. Dashed line: 
SBSE including the QP renormalization factors. Full line: result 
of the dynamical Bethe-Salpeter equation including dynamical QP 
and excitonic effects. ^Circles: experimental spectra (Ref. il7il for 
Cu and Ag and Ref. O for Si). 



The case of bulk silicon is also shown in Fig[21 bottom 
panel. The SBSE spectrum, calculated in the usual way 
by setting the Z factors to 1 hy hand {dotted line) is in 
good agreement with the experiment (circles). On the 
other hand, the spectrum including dynamical factors in 
the SBSE (dashed line) underestimates the experiment, 
as anticipated above, and the relative intensity of the 
two peaks is only poorly reproduced. The solution of 
the DBSE (solid line), instead, is in good agreement with 
the experiment. In contrast to the metallic case, how- 
ever, the DBSE kernel of silicon must contain second- 
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order contributions in order to reproduce correctly the 
experimental optical spectrum. The main effect of the 
first order kernel 11 {oj) is indeed to balance the reduc- 
tion of optical strengths due to self-energy renormaliza- 
tion factors, as suggested by Bechstedt et al. Q. How- 
ever, the renormalized QP weights also imply a reduction 
of the statically screened electron-hole of almost ~30%, 
which is the reason for the wrong relative intensities of 
the two peaks in the SBSE result. This shortcoming is 
fixed by the second-order diagrams, as discussed below. 
One of the two second order diagrams, pf^'") (t), is shown 
in Fig^ right diagram. In the other second-order di- 
agram (p(2'^) (t), not shown), the interaction lines are 
ordered according to t > i2 > ^4 > > ^3 > 0. Us- 
ing Eqs. H5I6|I we can isolate the external non-interacting 
polarization functions (as schematically shown in Fig^ 
by the crosses corresponding to the time points ti and 
^4). Consequently P^^.") (t) can be Fourier transformed 
to yield P^^^'^) {uj) = -P(o) (uj) n(2.'^) (w) p(o) (w) with 
n^^''') (uj) = 0('»)z^'*n('') (lu). Together with the contri- 
bution from diagram P^^-''), we obtain the total second- 
order kernel, 

n(2) (uj) = n (w) + 0'"^z"''n(*) (w) . (lo) 

s — a,b 

The solution of the DBSE with the kernel given by 
Eq. H1U|I is equivalent to summing all diagrams of the 
original BSE with up to two overlapping interaction 
lines. While the first-order term (11 (w)) partially re- 
stores the optical-strength intensities, the second-order 
correction (0('')Z'^''n'^*' (w)) reduces the screening of the 
electron-hole pair, thus enhancing their interaction. In 
the case of silicon this effect improves considerably the 
relative intensity of the two main peaks, as shown in 
Fig.|21(continuous line). Higher-order contributions can 
also be included in the DBSE kernel, although these 
terms contain either higher powers of or higher or- 
der energy derivatives of 11 (u) 19] that are negligible in 
the cases studied in this work. 

Thus, in the case of copper and silver the BSE must 
be solved using a time-dependent electron-hole interac- 
tion to obtain good agreement with the experiment. The 
DBSE result in the case of silicon is similar to the case 
when, in the SBSE, the QP renormalization factors are 
neglected j^- We understand this result as follows: the 
electron-hole pair, being a neutral excitation, is less ef- 
ficient than the electron and the hole alone in exciting 
virtual plasmons, which is the main process leading to 
QP renormalization. Only when dynamical effects are 
coherently included both in the self energy and in the 
electron-hole interaction, this (physically expected) re- 
sult emerges from the bundle of many-body equations. 
This confirms the SBSE results but not the separate ap- 
proximations involved therein. Even though the plasma 
frequency of silicon is at ^16 eV (far above the optical 



energy range), the DBSE kernel needs a second order con- 
tribution to balance the self-energy dynamical factors. 

In conclusion, we have shown that, in contrast to com- 
mon belief, dynamical excitonic effects in metals do exist 
and are crucial for reproducing the experimental optical 
absorption. We have demonstrated that dynamical exci- 
tonic and self-energy effects must be included together 
in the calculation of the response functions as they do 
not completely cancel each other. The good agreement 
with experiment obtained using a static BSE kernel for 
semiconductors is confirmed by the present, more general 
approach. 
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